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Abstract — This  paper  develops  five  data-aided  frequency  offset 
estimators  for  SOQPSK-TG  equipped  with  the  iNET  preamble 
and  operating  in  ISI  channels.  Four  of  the  five  estimators  exam¬ 
ined  here  are  generalizations  of  frequency  estimators  developed 
for  the  AWGN  channel.  Simulations,  performed  over  channel 
impulse  responses  measured  in  multipath  scenarios  encountered 
in  aeronautical  telemetry,  show  the  modified  L&R  estimator 
is  the  best  in  terms  of  tracking  range  and  achievable  mean- 
squared  error.  The  ESP  estimator  is  a  close  second  and  offers 
a  computational  complexity  advantage  in  that  it  only  requires 
the  computation  of  one  correlation  value  as  opposed  to  the  five 
correlation  values  required  by  the  L&R  estimator. 

I.  Introduction 

The  recent  introduction  of  the  iNET  (for  integrated  Net¬ 
work  Enhanced  Telemetry)  standard  [1]  provides  a  mode 
for  packetized  telemetry  downlinks.  Each  packet  comprises 
three  bit  fields:  a  preamble,  an  attached  sync  marker  (ASM), 
and  data  bits  (an  LDPC  codeword).  The  availability  of  a 
preamble  introduces  the  possibility  of  data-aided  synchro¬ 
nization  in  aeronautical  telemetry.  Compensating  for  carrier 
frequency  and  phase  offsets  and  timing  delay  offsets  is  re¬ 
quired  for  optimum  coherent  detection.  Multipath  interference, 
in  the  form  of  frequency  selective  fading,  complicates  the 
the  synchronization  problem.  Propagation  through  frequency 
selective  fading  introduces  inter-symbol  interference  (ISI)  into 
the  received  signal,  and  ISI  can  wreack  havoc  on  estimators 
designed  for  the  additive  white  Gaussian  noise  (AWGN) 
channel.  Eor  frequency  offset  estimation,  the  usual  approach 
is  to  use  a  preamble  with  some  periodic  structure.  With  proper 
preprocessing,  frequency  offset  estimators  designed  for  the 
AWGN  channel  may  then  be  applied,  such  as  the  maximum 
likelihood  estimator  formulated  by  Rife  and  Boorstyn  [2],  and 
the  reduced-complexity  estimators  described  by  Kay  [3],  Eitz 
[4],  Luise  &  Reggiannini  [5],  and  Mengali  &  Morelii  [6].  A 
nice  example  in  the  open  literature  is  from  harris  and  Dick 

[7]. 

In  this  paper,  we  examine  the  problem  of  estimating  the 
frequency  offset  in  an  ISI  channel  using  SOQPK-TG  with 
the  iNET  preamble.  We  exploit  the  periodic  structure  of  the 
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iNET  preamble  to  develop  a  low-complexity  frequency  offset 
estimator  with  good  performance. 


II.  Preliminaries 


The  bit  sequence  for  iNET  is  organized  as  outlined  in  the 
top  portion  of  Eigure  1 .  The  preamble  sequence  is  the  iNET 
standard  is  CD98hex  repeated  eight  times  [1,  p.  48].  This  bit 
pattern,  derived  from  the  analysis  presented  in  [8],  included 
the  repetition  to  enable  fast  acquisition.  Here,  the  periodic 
component  associated  with  the  repetition  is  leveraged  to  enable 
low-complexity  frequency  offset  estimation. 

The  transmitted  signal,  s{t),  is  SOQPSK-TG  whose  input 
bit  stream  is  summarized  in  the  top  portion  of  Eigure  1. 
SOQPSK-TG  is  a  partial  response  CPM  waveform  with  a 
constrained  ternary  alphabet.  The  details  are  spelled  out  in 
[9],  [10].  The  signal  propagates  through  a  frequency  selective 
channel  and  experiences  a  frequency  offset  and  the  addition  of 
additive  white  Gaussian  noise.  We  assume  the  received  signal 
is  filtered,  down-converted  to  complex  baseband,  and  sampled 
(not  necessarily  in  that  order)  using  standard  techniques.  The 
resulting  sequence  of  received  samples  is 


r(n) 


■  N2 

h{k)s{n  —  k) 

.k=^-Ni 


w{n) 


(1) 


where  h{n)  is  the  impulse  response  of  the  unknown  channel 
impulse  response  with  support  on  —  Vi  <  n  <  N2,  ojq 
rads/sample  is  the  frequency  offset  to  be  estimated,  and  w{n) 
is  a  proper  [11]  complex-valued  zero-mean  Gaussian  random 
process  with  auto-covariance  function 

^E{w{n)w*{n- k)}  =  a^Sik).  (2) 

Because  SOQPSK-TG  is  a  nonlinear  modulation,  the  estimator 
cannot  operate  on  the  symbols  in  the  same  way  it  does  for 
linear  modulation.  Consequently,  the  estimator  must  operate 
on  the  samples  of  SOQPSK-TG.  Let  i  be  the  index  in  r{n) 
corresponding  to  the  beginning  of  the  preamble  sequence  and 
let  Lp  be  the  number  of  samples  in  the  preamble  sequence. 
The  situation  is  illustrated  in  Eigure  1.  Because  each  sample 
depends  on  the  unknown  channel  impulse  response  h{n),  the 
standard  approaches  cannot  be  applied  directly.  However,  an 
estimator  that  leverages  the  periodic  properties  of  the  iNET 
preamble  can  be  applied. 
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Fig.  1.  A  graphical  illustration  of  the  iNET  data  packet  structure:  (top)  the  data,  preamble,  and  ASM  fields;  (bottom)  the  relationship  between  the  indexes 
of  the  received  samples  r(n),  the  signal  samples  s{n),  the  preamble  samples  p{n)  and  the  short  sequences  q{n)  that  constitute  the  preamble  sequence. 


Because  the  iNet  preamble  comprises  8  repetitions  of  a  16 
bit  sequence,  the  resulting  SOQPSK-TG  modulated  carrier  is 
also  periodic  over  the  time  interval  corresponding  to  the  oc¬ 
currence  of  the  preamble.  Furthermore  the  sampled  SOQPSK- 
TG  signal  is  periodic  over  the  interval  corresponding  to  the 
preamble  sequence.  The  situation  is  illustrated  in  Figure  1. 
Here,  s(n)  are  the  samples  of  the  complex  baseband  SOQPSK- 
TG  signal  and  p{n')  for  n'  =  0, 1, . . . ,  Lp  —  1  are  the  samples 
of  the  complex  baseband  SOQPSK-TG  signal  corresponding 
to  the  preamble  bits.  If  the  preamble  sequence  starts  at 
index  i,  then  s(z)  =  p(0),  s{i  -f  1)  =  p{l)  and  so  on  to 
s{i  +  Lp  —  1)  =  p{Lp  —  1).  Now  let  q{n)  be  the  samples  of 
the  complex  baseband  SOQPSK-TG  signal  corresponding  to 
the  16-bit  sequence  that  is  repeated  to  form  the  preamble,  and 
let  the  length  of  q{n)  be  Lg.  (At  an  equivalent  sample  rate  of 
2  samples/bit,  Lp  =  256  and  Lg  =  32.) 

For  n  =  i,i  +  1, . . .  ,i  +  Lp  —  1,  we  write  r{n)  =  r{i  +  b) 
for  6  =  0, 1, . . . ,  Lp  —  1  and  express  it  as  follows: 

■  N2 

r{i  +  b)  =  ^  h{k)s{i  +  b  —  k) 


=  h{k)p{b-k) 

.k=-Ni 

To  incorporate  the  periodic  nature  of  p(n)  into  the  equation, 
we  write 


-f  u;(z-f  5).  (4) 


(*+'’)+ u;(z  +  6) 


b  =  £Lg  +  m,  f  =  0, 1,...,7,  m  =  0, 1, . . . ,  L^  —  1.  (5) 
Here  we  see  that  I  indexes  the  block  and  m  indexes  the  sample 


within  the  block.  Substituting  gives 

r(i  -f  £Lg  -f  m)  = 

■  N2 

h{k)p{£Lq  +  m  - 

.k=-Ni 
■  N2 

=  h{k)q{m  —  k) 

.k=-Ni 

' - V - 

a{m) 

=  -f  w{i  +  (Lg  +  m).  (6) 

The  relationship  (6)  is  important.  It  shows  that  the  contribution 
of  the  channel  to  r{i  +  (Lg  +  m)  depends  only  on  to, 
the  position  inside  the  q{-)  block.  Consequently,  the  only 
difference  (other  than  the  additive  noise)  between  the  sample 
at  position  to  in  adjacent  blocks  is  the  phase  rotation  due  to  the 
frequency  offset.  ^  This  feature  is  captured  by  the  correlation 
function. 

Assuming  Ni  <  Lg  and  N2  <  Lg,  we  are  interested  in  the 
correlation  function  involving  the  middle  six  blocks: 

i-h7  Lq  —  1 

r{n)r*{n-S)  (7) 

^  n—i-\-Lq-\-S 

Because  the  samples  r(n)  are  periodic  with  period  Lg,  we 
are  interested  in  evaluating  this  function  at  integer  multiples 

^This  is  only  true  for  middle  six  blocks  (i.e.,  £  =  1, . . . ,  6)  because  the 
convolution  sum  in  the  first  block  {£  =  0)  includes  samples  from  the  data 
field  immediately  preceding  the  first  block  whereas  the  last  block  £  =  7 
includes  samples  from  the  ASM  field  immediately  following  the  last  block. 
Eliminating  the  first  and  last  blocks  is  sufficient  as  long  as  7Vi  <  Lq  and 
-/V2  <  Lg.  If  either  of  these  conditions  are  not  true,  then  (6)  is  true  for  fewer 
blocks. 


of  Lq.  For  S  =  dLq  (1  <  fi  <  5)  the  correlation 
Appendix) 

2 

R{dLq)  =  +  V 

where  is  given  by 

Lq  —  1 

^  \a{m)\^  (9) 

m— 0 

and  V  is  approximately  a  complex-valued  Gaussian  random 
variable  with  zero  mean  and  variance  (39). 

An  obvious  estimator  is  the  delay-and-multiply  (D&M) 
estimator  (cf  Eq.  (4.5.1)  of  [12]) 

Wo  =  ^arg|i?(Lq)|  .  (10) 

This  estimator  is  simple  and  exploits  the  noise  averaging  that 
occurs  over  the  5Lq  consecutive  samples  of  r(n).  But  this 
estimator  only  uses  R{Lq).  It  would  seem  there  might  be 
performance  advantages  to  incorporating  the  information  in 
R{2Lq),  R{‘iLq),  R{ALq),  and  R{5Lq).  To  see  how  this  might 
be  done,  we  first  look  to  how  this  same  concept  is  applied  in 
the  AWGN  case.  Next  we  modify  the  developments  to  create 
frequency  estimators  for  present  case. 

III.  Data-Aided  Frequency  Estimators  eor  AWGN 

The  standard  model  for  data-aided  frequency  estimators  in 
the  additive  white  Gaussian  noise  environment  is  the  sinusoid 
in  noise. 

x{n)  =  +  w{n)  (11) 

for  n  =  0, . . . ,  N  —  1.  The  sinusoid  in  noise  model  is 
usually  obtained  from  the  sampled  matched  filter  outputs  at 
the  receiver  by  multiplying  each  matched  hlter  output  by  the 
conjugate  of  the  corresponding  pilot  symbol.  The  maximum 
likelihood  (ML)  estimator  is  [2] 


r 

A^-1 

cDo  =  argmax  < 

x{n)e-^‘^°^ 

[ 

( 

1 

n—0 

1 

In  words,  the  ML  estimate  is  the  one  that  maximizes  the 
periodogram  of  the  sequence  x{n).  Because  there  is  no  closed 
form  solution,  this  approach  involves  a  search.  For  computa¬ 
tional  complexity  reasons,  a  suboptimal  solution  that  does  not 
require  a  search  is  desirable.  The  most  common  data-aided 
estimators  operate  on  the  x{n)  or  the  autocorrelation  function 
of  the  x{n). 

Fitz  Estimator:  The  Fitz  estimator  is  [4] 

2  ^ 

=  M(M  +  1)  S  H  {«-(’»)})  03) 

where  R^im)  is  the  unnormalized  autocorrelation  function 

N-l 

Rxini)  =  x{n)x*{n  —  m).  (14) 


is  (see  the 
(8) 


L&R  Estimator:  The  Luise  and  Reggiannini  (L&R)  estimator 
is  [5] 


M 


Wo  = 


M  +  1 


arg  E  Rxijn) 


(15) 


where  Rx{m)  is  the  normalized  autocorrelation  function 


N-l 


^  ^ 

(tit)  =  -  /  x{n)x*{n  —  m).  (16) 

N  —  m 


M&M  Estimator:  The  Mengali  and  Morelli  (M&M)  estimator 
is  [6] 


Here,  1  is  the  M  x  1  all-ones  vector  and 

0  =[</>(!)  m  ■■■  m)] 


(18) 


where 


(j){m)  =  arg  {i?^(m)}  -  arg  {R^{m-  1)}  (19) 


is  the  phase  difference  between  successive  autocorrelation 
functions  (16).  (The  phase  differences  are  used  to  overcome 
the  phase  wrapping  problem  with  arg{-}.)  C  is  the  M  x  M 
covariance  matrix  of  the  noise  samples  associated  with  the 
autocorrelation  phase  differences. 

Kay  Estimator:  The  Kay  estimator  is  [3] 


Wo 


l^C-^A 


(20) 


where  1  is  the  (A^—  1)  x  1  all-ones  vector,  A  is  the  (A^—  1)  x  1 
vector 


A=[A(1)  A(2)  •••  A(A-l)]^  (21) 


where 

A(n)  =  arg  {a;(n)a;*(n  —  1)}  (22) 


is  the  phase  difference  between  adjacent  x{n),  and  C  is  the 
{N  —  1)  X  {N  —  1)  covariance  vector  of  the  noise  samples 
accompanying  the  phase  differences. 


IV.  Extensions  to  iNET  Preamble-Based 
Frequency  Estimators  eor  ISI  Channels 

In  the  previous  section,  we  saw  that  the  Fitz,  L&R,  and 
M&M  estimators  used  information  in  the  autocorrelation 
function  (both  unnormalized  and  normalized)  of  the  received 
data  samples.  The  modification  of  the  Fitz,  L&R,  and  M&M 
estimators  more  complicated  than  simply  substituting  R{dLq) 
[see  (7)]  for  Rx{m).  This  is  because  the  sample  spacing 
between  the  R{dLq)  is  Lq  samples  instead  of  one  sample  for 
Rx{m).  This  not  only  impacts  the  noise  correlations  but  also 
modihes  the  approximations  applied  in  the  intermediate  steps 
in  the  derivations.  Because  space  limitations  do  not  permit  a 
detailed  accounting  of  all  these  steps,  only  the  hnal  results  are 
listed  here. 


Modified  Fitz  Estimator:  The  modified  Fitz  estimator  is 


10 


1  ^ 

Wo  =  ^  marg  (23) 

^  m— 1 

where  R{mLq)  =  ((6  —  m)Lq)R{rnLq)  where  R{mLq)  is 
given  by  (7)  for  m  =  1,  2,  3, 4,  5. 

Modified  L&R  Estimator:  The  modified  L&R  estimator  is 

f  5  'i 


Wo  =  ^  arg  ^  R{mLq 
^  I  771=1 


(24) 


R{mLq)  is  given  by  (7)  for  m  =  1,  2,  3, 4,  5. 

Modified  M&M  Estimator:  The  modified  M&M  estimator  is 
1 


Wo  = 


where 


35L„ 


b(l){Lq)  +  8<^(2L,)  +  9<^(3L,) 

+  8</.(4L,)  +  5(^(5L,) 


(25) 


(j){mLq)  =  arg  |i?(TOLq)|  —  arg  [R{{m-l)Lq)}  (26) 

for  m  =  1,  2,  3, 4,  5. 

Modified  Kay  Estimator:  The  modified  Kay  estimator  is 


Wo  = 


Lq  —  1 

E 


a{m)\ 


35Lq  ^ 

y  771  =  0 


-  ^5Am,2  +  8Am,3 

+  8Am,5  +  5Am,6 


(27) 


where 

=  arg  |r(i  +  iLq  +  m)r*{i  +  {£  -  l)Lq  +  to)|  . 


The  challenge  with  this  form  of  the  estimator  is  that  the  a{m) 
are  unknown.  The  Q!(m)  are  used  to  weight  each  block  of 
phase  differences  (and  this  weighting  is  in  proportion  to  the 
relative  channel  output  power).  A  suboptimal  approach  is  to 
set  the  weights  to  be  equal: 


|q!(to)P  ^  1 

q;2  Lq 


(28) 


The  result  is 


1  «  / 

Wo  =  nrj^2  E  +  8Am_3 

m=0 

+  9Am,4  +  8Am,5  +  SA^.O^  •  (29) 
V.  Simulation  Results 

The  estimators  were  applied  to  three  representative  channel 
impulse  responses  capture  during  channel  sounding  exper¬ 
iments  at  Edwards  AFB,  California  [13].  These  channels 
are  illustrated  in  Figure  2.  See  the  figure  caption  for  more 
details.  In  these  simulations,  the  iNET  preamble  described  in 
Section  II  was  used  with  SOQPSK-TG  operating  at  10.3125 
Mibts/s,  the  “over-the-air  rate”  corresponding  to  a  payload  bit 
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frequency  (MHz) 


Fig.  2.  The  example  channels  from  channel  sounding  experiments  at  Edwards 
AFB:  (top)  a  length-9  channel  from  the  flight  line;  (middle)  a  length-19 
channel  from  take-off;  (bottom)  a  length-5  channel  from  low-elevation  angle 
“up  and  away”  flight.  In  each  plot,  the  thick  line  is  the  channel  frequency 
response  and  the  thin  line  is  the  power  spectral  density  of  SOQPSK-TG 
operating  at  10.3125  Mbits/s. 


rate  of  10  Mbits/s.  The  channel  and  SOQPSK-TG  preamble 
signal  operated  at  a  sample  rate  equivalent  to  2  samples/bit. 

As  the  performance  measures,  we  use  the  means  and 
variances  of  the  Penins  (ESP),  Fitz,  L&R,  M&M,  and  Kay 
estimators  over  the  three  channels  displayed  in  Figure  2.  The 
simulation  results  are  summarized  in  Figures  3-8.  Figures  3  - 
5  show  the  simulated  mean  of  the  five  estimators  over  the  three 
channels,  respectively.  All  five  perform  reasonably  well  up  to 
±60  kHz.  Here,  the  Fitz  and  M&M  estimators  stop  tracking. 
The  ESP,  L&R,  and  Kay  estimators  track  up  to  ±100  kHz 
(roughly  1%  of  the  bit  rate),  except  the  Kay  estimator  shows 
some  bias.  This  bias  channel  dependent  and  are  traced  back 
to  inaccuracies  in  the  approximation  (28). 

The  simulated  variances  of  the  estimators  are  summarized 
in  Figures  6-8  for  the  three  channels,  respectively.  These 
results  show  that  the  Fitz,  L&R,  and  M&M  estimators  achieve 
the  same  mean-squared  error  performance.  The  ESP  estimator 
is  a  few  dB  worse  than  these  three  whereas  the  Kay  estimator 
has  the  worst  mean-squared  error  performance.  This  is  not  a 
surprising  result  for  the  Kay  estimator  given  the  fact  that  the 
approximation  (28)  is  required. 


VI.  Conclusions 

On  the  whole,  the  modified  L&R  estimator  is  the  best  in 
terms  of  tracking  range  and  achievable  mean-squared  error. 
The  ESP  estimator  is  a  close  second  and  offers  a  computational 
complexity  advantage  in  that  it  only  requires  the  computation 
of  one  correlation  value  as  opposed  to  the  five  correlation 
values  required  by  the  L&R  estimator. 
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Appendix 

Here,  the  derivation  of  (8)  is  outlined.  Starting  with  the 
definition  (7)  and  using  6  =  dLg  gives 

i-\-7Lq  —  l 

R{dLg)  =  .g  _  XI  r{n)r*{n-5)  (30) 

^  n—i-\-Lq-\-5 

Using  the  indexes  (5),  R{dLq)  may  be  re-written  as 


R{dLg)  = 

6  Lq-1 

E  E  r  (i  +  ILq  +  'm)r*{i  +  {i—d)Lq+m). 

(0  - 


(31) 


Substituting  the  relationship  for  r{n)  given  in  (6)  and  rear¬ 
ranging  gives 


R{dLq)  = 
1 

(6  -  d)Lq 


6  Lq  —  1 

E  E 


2^3dL,uo 


(6  —  d)LQ 
6  Lq  —  1 

E  E[  Vi  {£,  m)  +  V2  {£,  m)  +  (f ,  m) 


m—Q 


The  first  term  in  (32)  may  be  simplified  as 


1 


6  Lq  —  ^ 


X  X 


2  jdLqUJo  _ 


where 


(6  —  d)Lq 

''  ’  1  t=d+l  m=0 


Lq-1 


=  X  l«Mr 


(32) 


(33) 


(34) 


m—0 


The  second  term  in  (32)  involves  the  double  summation  of 
three  noise  terms.  The  first  of  these  noise  terms  is 

vi{£,  m)  =  w* {i  +  {£  —  d)Lq  +  m)  (35) 


which  is  a  complex-valued  Gaussian  random  variable  with 
variance  \a(m)\‘^a^  [see  (2)].  Similarly, 

V2{£,m)  =  a* +  £Lq  +  m)  (36) 

is  also  a  complex-valued  Gaussian  random  variable  with 
variance  \a{m)\‘^a'^.  Furthermore,  given  the  assumption  that 
the  w{n)  are  a  sequence  of  uncorrelated  proper  complex¬ 
valued  Gaussian  random  variables,  vi  (£,  m)  and  V2  {£,  m)  are 
uncorrelated  [11].  The  last  of  the  noise  terms  in  (32)  is 


v^{£,  m)  =  w{i  -\-  £Lq  +  m)w*  {i  +  {£  —  d)Lq  +  m).  (37) 

Here  we  make  the  assumption^  that  for  sufficiently  high  signal- 
to-noise  ratios  (i.e.,  sufficiently  small  ctJ,),  this  product  term 
is  much  smaller  than  the  first  two  terms  with  high  probability. 
Consequently,  the  product  term  may  be  neglected  and  we  are 
left  with 


X  X  fL^)  +  V2{£,m) 


(38) 


which  is  the  sum  of  uncorrelated  complex-valued  Gaus¬ 
sian  random  variables  each  with  zero  mean  and  variance 
|q;(to)P(tJ.  Consequently,  v  is  complex-valued  Gaussian  ran¬ 
dom  variable  with  zero  mean  and  variance 

'  (39) 


a2  = 


LI 


^This  is  a  standard  assumption  in  the  frequency  estimation  literature.  See 
[3],  [4],  [5],  [6]. 
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Fig.  3.  Estimator  means  for  the  first  channel  of  Figure  2.  For  these 

simulations  E^j/Nq  =  15  dB  Estimator  variances  for  the  first  channel  of  Figure  2. 
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Fig.  4.  Estimator  means  for  the  second  channel  of  Figure  2.  For  these 

simulations  Eij/Nq  =  15  dB  Estimator  means  for  the  second  channel  of  Figure  2. 
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Fig.  5.  Estimator  means  for  the  third  channel  of  Figure  2.  For  these 

simulations  Ei^/Nq  =  15  dB  Estimator  means  for  the  third  channel  of  Figure  2. 


